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1 SPINLESS SALPETER EQUATION 

The conceptually simplest bound-state wave equation incorporating to some extent relativistic effects 
is the "spinless Salpeter equation." The origin of this equation of motion is by no means some mystery. 
Rather, this spinless Salpeter equation has to be regarded as a well-defined standard approximation to 
the Bethe-Salpeter formalism. 

1.1 Bethe— Salpeter (BS) Formalism 

The appropriate framework for the description of bound states within some relativistic quantum field 
theory is the Bethe-Salpeter formalism. Consider the bound state |M) (with momentum K and energy 
E) of, say, (spin-i) fermion and antifermion (with masses Mi, M 2 and momenta P lt P 2 , respectively). 
Within the BS formalism this bound state is represented by its Bethe-Salpeter amplitude ^ which, in 
momentum space, is defined as Fourier transform of the time-ordered product of the field operators of 
the two bound-state constituents between the vacuum |0) and the bound state, after factorizing off the 
centcr-of-momcntum motion: 



*(P) = exp(i^XcM) / d 4 X exp(iPX) (0|T(V>i(Xi)^(X 2 ))|M(lO} 



with the usual kincmatical notions of total momentum K = P1+P2, relative momentum P, center-of- 
momentum coordinate Xcm, and relative coordinate X = X\ — X 2 of the bound-state constituents. 
The BS amplitude W satisfies the Bethe-Salpeter equation pi 

S^ 1 (P 1 )^(P)S 2 1 (-P 2 ) = -^ J d 4 Q/C(P,Q)*(Q) , 



which involves two dynamical ingredients, namely, 

• on the one hand, the full fermion propagator of particle i (i = 1,2), denoted by S(P, M) for any 
particle of momentum P and mass M, which, however, is usually approximated by its free form 

S^(P,M) = -i( ltl P^-M) , 

with M interpreted as some effective ( "constituent" ) mass and the propagator So understood as 
an effective one, and, 

• on the other hand, the BS interaction kernel K.(P, Q), which is defined (only perturbatively!) as 
sum of all two-particle ("BS-") irreducible Feynman diagrams for two-particle into two-particle 
scattering (exhibiting, consequently, a certain amount of gauge dependence), and which is given, 
for instance, in lowest order quantum electrodynamics from one-photon exchange (in Feynman 
gauge) by 

By construction, the BS equation is formally exact. Its actual application, however, faces several, 
rather well-known problems]^] 

• There is no means to compute the BS kernel beyond the very tight limits of perturbation theory. 

• Even with the BS kernel at one's disposal, it is, in general, not possible to find the exact solutions 
of the BS equation (except for rather few special cases, like the famous Wick-Cutkosky model 
describing the interaction of two scalar particles by exchange of some massless scalar particle). 

• In non-Abelian gauge theories like, e.g., quantum chromo dynamics, since the Dyson— Schwinger 
equations connect the propagators and the interaction kernel, the simultaneous assumption of 
free propagators and a confining BS kernel is obviously inconsistent. 

The most straightforward way out is the reduction of the Bethe-Salpeter equation by a series of 
approximations: 

1 The appearance of states with negative or vanishing norm among the obtained solutions as well as 
questions of interpretation, for instance, cause additional troubles when dealing with the BS equation. 
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1. Eliminate any dependence on timelike variables. This (innocent) enterprise involves two steps: 

• Adhere to the static approximation to the BS kernel (which has to be justified a posteriori) , 
that is, assume that the BS kernel depends only on the relative three- momenta Q and P 
of initial and final state, respectively: 

K{P,Q) =/C(P,Q) . 

(This amounts to ignoring all retardation effects by assuming instantaneous interactions.) 

• Define the equal-time wave function $, sometimes also called the "Salpeter amplitude," (in 
momentum space) by integrating the BS amplitude \&(P) over the zero component of the 
four-momentum P: 



$(P)= / dP *(P,P ) 



This leads to the Salpeter equation 



<D(P) 



d 3 Q 
(2tt) 3 



A+ 7o /C(P, Q) $(Q) 7 o (A+) c A^ 7 o £(P, Q) *(Q) 7o (Aj 



E 



P? 



Ml 



PI 



Ml 



E- 



P? 



M 2 



PI 



with the energy projection operators for positive or negative energy of particle i (i = 1, 2) given, 
as usual, by 



Mf± l0 ( 1 -P l +M l ) 



Mr 



and 



Af 



where C is the charge conjugation matrix. The Salpeter equation is the equation of motion for 
the equal-time wave function $(P), with full relativistic kinematics but in static approximation 
for the interaction kernel JC. 



2. Neglect the second term on the right-hand side of the Salpeter equation on the basis of the — at 
least for "heavy" bound-state constituents reasonable — assumption that the denominator of the 
first term is much smaller than the denominator of the second term: 



P? 



Ml 



Ml<^E 



P? 



Ml 



PI 



Ml 



This leads to the reduced Salpeter equation 



E 



P? 



Ml 



Pi 



Ml $(P) 



d 3 Q 
(2tt) 3 



A+ 7o /C(P,Q)$(Q)7o(A+; 



3. Neglect any reference to the spin degrees of freedom of the involved bound-state constituents. 

4. Assume that the BS kernel JC(P, Q) depends only on the difference of the relative momenta P 
and Q, which means that the BS kernel is of convolution type (as is, in fact, frequently the case): 

/C(P,Q) =/C(P-Q) . 



5. Restrict the whole formalism exclusively to positive-energy solutions, which will be denoted by 

After applying all these simplifying assumptions and approximations to the BS equation, one finally 
ends up with the spinless Salpeter equation 



Pf + Ml 



Pi 



Ml 



V(X) 



ip = Eip 



involving an interaction potential, V(X), arising as the Fourier transform of the BS kernel JC(P — Q). 
In particular, in the center-of-momcntum frame of the two bound-state constituents (i.e., for K = 0), 
this equation reads 

H\*l>) = E\i>) , 

with the Hamiltonian 

H = 



P 2 + Ml 



P 2 + Ml + V(X) 
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1.2 Equal-Mass Case 

For the special case of equal masses of the two involved bound-state constituents, i.e., when assuming 

M x = M 2 = M , 

it is possible to recast the two-particle spinless Salpeter equation into the equivalent one-particle form. 
In order to see this, consider the semirelativistic Hamiltonian H for the two-particle spinless Salpeter 
equation with an interaction described by a central potential t^(|X|) of, for instance, power-law form: 

H = 2 VP 2 + M 2 + fc ™ Rn > R = l X l • 

One may always perform a scale transformation of the phase-space variables X, P by some arbitrary 
scale factor A, 

p = AP, x=y, 

since any transformation of this kind (necessarily) preserves the fundamental commutation relations: 

[x,p] = [X,P] . 

Now, 

• fix the scale parameter A to the value A = 2, which implies for the semirelativistic Hamiltonian 

H = V / 4P 2 + 4M 2 + k n R n 

= v / P 2 +4Af 2 + ^fc„2"r" , 

• identify the mass and coupling-strength parameters for one- and two-particle form according to 

m = 2 M , 
a n = 2" k n , neZ, 
and arrive — presumably without great surprise — at the equivalent one-particle Hamiltonian 

H = v/p 2 + m 2 |^a„r", r = |x| . 

nGZ 

1.3 (One-Particle) Spinless Salpeter Equation 

In view of the above observation, it is sufficient to confine the present discussion to the consideration of 
some self-adjoint Hamiltonian H of the form 

H = T + V, (1) 

where T denotes the "square-root" operator of the relativistic expression for the free (kinetic) energy 
of some particle of mass m and momentum p, 

T = T(p) = VP 2 + m 2 , 

and V = V(x) represents some arbitrary coordinate-dependent, static interaction-potential operator. 
The spinless Salpeter equation is nothing else but the eigenvalue equation for the Hamiltonian H, 

H\ X k)=E k \ Xk ) , k = 0,1,2,... , 

for the (complete set of) Hilbert-space eigenvectors \xk) of H corresponding to the energy eigenvalues 

Ek ^ (Xk\H\xk) 

(Xk\Xk) 

As such, it represents the simplest relativistic generalization of the Schrodingcr equation of standard 
nonrelativistic quantum theory. N.B. The semirelativistic Hamiltonian H is a nonlocal operator, i.e., 

• either the relativistic kinetic-energy operator T in configuration space 

• or, in general, the interaction-potential operator V in momentum space 

is nonlocal. Because of this nonlocality it is somewhat difficult to obtain rigorous analytic statements 
on the solutions of this equation of motion. 
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2 RELATIVISTIC COULOMB PROBLEM 

Of particular importance in many realms of physics is the (spherically symmetric) Coulomb potential, 
with interaction strength parametrized by some dimensionless coupling ("fine structure") constant a: 

V(x) = V c (r) = -- , r=|x|, a >0 . (2) 
r 

2.1 Spinless Relativistic Coulomb Problem 

Since we will use the Coulomb potential in order to illustrate the general statements derived below, let 
us briefly summarize, in roughly chronological order, the most important knowledge gained until now: 

• By examination of the spectral properties of the semirelativistic Hamiltonian H with a Coulomb 
interaction potential Vc (r) , one is able to prove || for this operator H 

— its essential self-adjointness (which means that the closure of H is self-adjoint) for a < — ; 

— the existence of its Friedrichs extension up to a critical value a c of the coupling constant a: 

2 

a c = - ; 

TT 

— a strict lower bound on the ground-state energy Eq (that is, on the operator H), given by 



TT a \ ^ 2 

Eq > m \l 1 — ( —^-j for a < — 



Loosely speaking, the Coulombic Hamiltonian H may be regarded as reasonable operator up to 
the critical coupling constant a c . 

For part of the allowed range of the coupling constant a, an improved lower bound can be found 




Eq > m 

• The analytic solution for the corresponding wave function ip has been constructed |5j ; inspecting 
its behaviour near the origin, one observes for, e.g., vanishing orbital angular momentum [I = 0) 
that the configuration-space wave function tp(x) diverges for small values of r = |x| exactly like 

hm wx) oc r with v = — 1 H 

• Given the relative uselessness of Temple's inequality, lower bounds on the spectrum of H may be 
found (in a numerical manner) with the help of the generalized "local-energy" theorem 0: 

— Assume 

1. that the Fourier transform V(p) of the interaction potential V(x) is strictly negative, 
except at infinity, as is (certainly) the case for the (attractive) Coulomb potential (g), 

2. that the spectrum of the operator H is discrete, and 

3. that the ground state of the Hamiltonian H exists. 

— Define the "local energy" 



£(P) = T(P) 



J d 3 qV(p-q) <Mq) 



<P(P) 

where </>(p) denotes some suitably chosen, positive trial function, 

<Mp) > o . 

Then the lowest- lying eigenvalue Eq of the Hamiltonian H is bounded from below and above by0 



inf£(p) < E Q < sup£(p) 
p P 



2 The lower bound even holds if the assumption on the discreteness of the spectrum is not fulfilled. 
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Table 1: Numerically (!) obtained Q lower and upper bounds on the ground-state energy E (a) 



a — 

m 



Lower Bound Upper Bound 



0.0155522 


0.9998785 


0.9998791 


0.1425460 


0.989458 


0.989613 


0.2599358 


0.96309 


0.96364 


0.3566678 


0.92578 


0.92673 


0.4359255 


0.88013 


0.88139 


0.5 


0.82758 


0.82910 


0.5505528 


0.76908 


0.77075 


0.5887832 


0.7050 


0.7069 


0.6155367 


0.6359 


0.6379 


0.6313752 


0.5616 


0.5637 


0.6366198 


0.4825 


0.4843 



With the help of this local-energy theorem, the ground-state energy, considered as a function Eq{q) of 
the coupling strength a, has been restricted to some remarkably narrow band (Table 0). In particular, 
at the critical coupling constant a c one finds 

-En (a = a r ) 

0.4825 < — < 0.4842910 for a = a c . 

to 

From this one learns that the ground-state energy E at the critical coupling constant a c is definitely 
nonvanishing (which is evidently not clear from the lower bounds quoted above). 

2.2 Power Series Expansions 

The systematic series expansions of the energy eigenvalues E(n T , £) in powers of the coupling constant 
a (up to and including the order 0(a 7 )) for states with radial quantum number n r and orbital angular 
momentum I have been worked out in different ways [0, || . The disadvantage of these power-series 
expansions is that they are only significant for a region of rather small values of the coupling constant 
a. These expansions read, for instance, 

• for n T = 0, I = (i.e., for the ground state): 



E(n r = 0,£ = 0) _ a 2 5 a 



4 



= 1 1 h a b In 



to 2 8 3vr 

7 2 33 

^2 ~ 16 



n 



C(3) , 



) a 6 + 0{a 7 ) , 



where is the Riemann zeta function 

oo 

C(z) = Y / k- z , Kz>l, 
fc=i 

which coincides up to order 0(a 4 ) with the improved lower bound on the ground-state energy; 

• for n r = 1, I = 0: 

E(n T = 1,£ = 0) a 2 13 a 4 a 5 a 6 lna 
_LJ ! L = i 1 1 

to 8 128 3tt 8 

8tt 2 SV ; 4tt 2 1024 8 J K ' ' 

• for n T = 0, t = 1, where, remarkably, there are no contributions of the order 0(a 5 ) or 0(a 6 In a): 

E(n T = 0,£=1) a 2 7a 4 727a 6 ^. 7 , 
to = 1 "T- 384 + 82944 + ° (a } ' 

In any case, even for the Coulomb potential (Q) the eigenvalues of the semirelativistic Hamiltonian (|l|) 
are not known exactly! 
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3 EXACT ANALYTIC UPPER BOUNDS ON ENERGY 
LEVELS 

In view of the (admittedly rather unsatisfactory) state-of-the-art mentioned above, let us try to find, 
at least, exact upper bounds on the eigenvalues of the Hamiltonian ([j]) by entirely analytic methods. 
Clearly, the derivation of upper bounds on the eigenvalues of some operator only makes sense if this 
operator is bounded from below. Accordingly, let us assume that the (otherwise arbitrary) interaction 
potential V is such that the Hamiltonian H is bounded from below. For the Coulomb potential (0), 
the required semi-boundedness of the spectrum of the Hamiltonian H has been rigorously proven |3j . 



3.1 Minimum— Maximum Principle 



Beyond doubt, the theoretical foundation of (as well as the primary tool for) the derivation of rigorous 
upper bounds on the eigenvalues of some self-adjoint operator is the well-known "min-max principle." 
There are several equivalent formulations of this theorem. For practical purposes, the most convenient 
one is the following: 

• Let H be a self-adjoint operator bounded from below. 

• Let Ek, k = 0, 1, 2, . . ., denote the eigenvalues of H, 

H\ Xk )=E k \xk) , fc = 0,1,2,... , 

ordered according to 

E < E x < E 2 < . . . . 

• Let Dd be some d-dimensional subspace of the domain of H . 

Then the fcth eigenvalue Ek (counting multiplicity) of H satisfies the inequality 



E k < sup 

WeD k+1 



(jW) 
(V#> 



for k = 0,1,2,. 



3.2 Operator Inequalities 

Needless to say, our goal must be to replace the problematic kinetic-energy square-root operator in the 
semirelativistic Hamiltonian H by some more tractable operator. One (obvious) way to achieve this is 
to use the min-max principle in order to compare eigenvalues of operators: 

• Assume the validity of a generic operator inequality of the form 

H < O . 



Then 



Ei- 



_ (Xk\H\ Xk ) 

(Xk\Xk) 

{MHU) „ , . . . , . 

< sup — - — j — - — (by the mm-max principle) 



(V#> 



< sup — — j — — (by the operator inequality) . 

• Assume that the (fc+ l)-dimensional subspace -Dfc+i in these inequalities is spanned by the first 
k+1 eigenvectors of the operator 0, that is, by precisely those eigenvectors of O that correspond 
to the first k + 1 eigenvalues Eq, Ex, . . . , of O if all eigenvalues of O are ordered according to 

Eq ^ Ex ^ Ei 5^ • • • ■ 

Then 

sup M = £fc . 

Men k+1 (V#> 

Consequently, every eigenvalue Ek of H is bounded from above by a corresponding eigenvalue Ek of O: 

Ek < Ek ■ 
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3.3 The "Schrodinger" Bound |g| 

The simplest of all upper bounds on the eigenvalues of the semirelativistic Hamiltonian H of Eq. ([!]) is 
the one involving the corresponding "Schrodinger" Hamiltonian. It may be found by taking advantage 
of the positivity of the square of the (since T is self-adjoint, obviously self-adjoint) operator T — m: 

< {T-mf 

= T 2 + m 2 -2mT 
= p 2 + 2m 2 -2mT. 



Assuming m to be positive, this may be converted into an operator inequality for the kinetic energy T, 

2 m ' 



1 < to 



which, in turn, entails an operator inequality for the (generic, since it involves an arbitrary interaction 
potential V) Hamiltonian H in the spinless Salpeter equation |5): 

H <H S , 

with the "Schrodinger" Hamiltonian 

P 2 

2 TO 

For the Coulomb potential (0), the energy eigenvalues corresponding to the above Schrodinger 
Hamiltonian read 

a 2 



Es, n = m 1 

\ 2 n z 

where the total quantum number n is given in terms of both radial and orbital angular-momentum 
quantum numbers n r and £, respectively, by 

n = n T + £+l, n r = 0,1, 2, ... , £ = 0,1,2,.... 
3.4 A "Squared" Bound 

In order to improve the above bound, one might be tempted to consider the square of the Hamiltonian 
H: 

Q = H 2 =T 2 + V 2 + TV + VT. 



The eigenvalue equation for this squared Hamiltonian Q will, of course, be solved by exactly the same 
set of eigenvectors |xfe) as the one for the original Hamiltonian H with, however, the squares E 2 of the 
corresponding energy eigenvalues Ek of H as the eigenvalues of Q: 

Q\Xk) = El\ Xk ) , k = 0,1,2,.... 

The positivity of the square of the self-adjoint operator T — m — V, 

< (T-m- V) 2 

= T 2 + m 2 + V 2 - 2mT+2mV -TV - VT , 

and the (because of the positivity of the operator p 2 obviously valid) relation 

< to < T 

yield some operator inequality for the anticommutator TV + VT of kinetic energy T and interaction 
potential V showing up in the squared Hamiltonian Q: 

TV + VT < T 2 + m 2 + V 2 -2mT+2mV 
= p 2 + 2 to 2 + V 2 - 2 to T + 2 to V 
< p 2 + V 2 + 2mV. 
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Insertion of this intermediate result implies an operator inequality for the squared Hamiltonian Q j^] : 

Q<R, 

with the operator 

R = 2 p 2 + m 2 + 2 V 2 + 2 m V . 

According to the min-max principle, the squares of the energy eigenvalues E k of the spinless Salpeter 
equation are therefore bounded from above by the corresponding eigenvalues £r^ of the operator R: 

E k < £rm 



Ek < y/£R,k ■ 

Only for the Coulomb potential (0), the above operator R is of exactly the same structure as the 
"Schrodinger" Hamiltonian H$, with, however, some "effective" orbital angular momentum quantum 
number L to be determined from the relation 

L(L + 1) = £(£+!) + a 2 , 

with the result 

Jl + 4 [I Ti + 1) + a 2 } - 1 
L= V L^_J J , £ = 0,1,2,.... 

The eigenvalues £ r of the operator R may thus be found by replacing in the Coulomb eigenvalues Es jU 
the angular momentum quantum number £ by the effective angular momentum quantum number L: 

„2 

2 



£r,n = m \ 1 



2N 2 . 

with the "effective" total quantum number N defined by 

N = n r +L + 1 , n T = 0,1,2,... . 
Unfortunately, these "squared" bounds lie above and are thus worse than the "Schrodinger" bounds. 

3.5 A Straightforward Generalization [|9[] 

The "Schrodinger" bounds may be improved by generalizing a little bit the line of argument leading to 
their derivation. The positivity of the square of the (obviously self-adjoint) operator T — ji, where /i is 
an arbitrary real parameter (with the dimension of mass) , 

< (T-fi) 2 

= T 2 + [i 2 - 2 n T 
= p 2 +TO 2 +/i 2 -2^T, 
implies a set of operator inequalities for the kinetic energy T|^| 

r< p 2 +m 2 + M 2 

This immediately translates into a set of operator inequalities for the semirelativistic Hamiltonian H 
§: 

H<H s (n) for all fi > , 
with the Schrodinger-like Hamiltonian 

2 \i 

Invoking again the min-max principle, any energy eigenvalue of H is thus bounded from above by 
the corresponding eigenvalue -Es./c(m) °f the Schrodinger-like Hamiltonian Hg(/j,) for given positive fi, 

E k <E Stk (n) forall^>0, 
and — as a (trivial) consequence — also by the minimum of all these Schrodinger-like upper bounds H : 



E k < minS s ,fe(A*) 

M>0 



3 See also Ref. g. 
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For the Coulomb potential (En, it is trivial to write down the corresponding energy eigenvalues: 



2/j, 



m 2 + /i 2 I 1 



with the total quantum number 



n = Th + 1 + 1 • 
Minimizing _Es,n(A*) with respect to the parameter yields 



minion (m) = m\ I 

M>0 



for all a < a c 



(3) 



These bounds 



• hold for all values a < a c of the Coulomb coupling constant a and arbitrary levels of excitation, 
and 

• improve definitely the Schrodinger bounds, for any value of the total (or "principal") quantum 
number n: 

min£ , s ! „(^) < E Sn for a ^ . 

(For fj, = m one necessarily recovers the Schrodinger case.) The comparison of the upper bounds of 
Eq. (|^) with their numerically Q or by a perturbative expansion |^] obtained counterparts shows that 

1. for the ground state, i.e., n r = I — 0, hence n = 1 (Table ||), the relative error of our analytical 
(operator) bound (in the considered range of the coupling constant a) is less than 4.5 %, and 

Table 2: Numerical |3| vs. analytical upper bounds on the ground-state energy Eq 





Upper Bound 
Numerical 


E 
on — 

m 

Analytical 


0.0155522 


0.9998791 


0.9998791 


0.1425460 


0.989613 


0.989788 


0.2599358 


0.96364 


0.96563 


0.3566678 


0.92673 


0.93423 


0.4359255 


0.88139 


0.8998 


0.5 


0.82910 


0.86603 



2. for the level n r = 0, £ — 1 and thus n = 2 (Table ||) the relative error of our analytical (operator) 
bound (in the considered range of the coupling constant a) is less than 0.1 %. 

Table 3: Perturbatively computed Iq] energy eigenvalues vs. analytical upper bounds for n r = 0, 1 = 1 



E(n T =0,^ = 1) 

a 

m 

Perturbation Theory Analytical Bound 



0.0155522 
0.1425460 
0.2599358 
0.3566678 
0.4359255 
0.5 



0.999969765 

0.997452 

0.99147 

0.9838 

0.975 

0.967 



0.999969766 

0.997457 

0.99152 

0.9840 

0.976 

0.9682 
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3.6 Rayleigh Ritz Variational Technique 

An immediate consequence of the min max principle is the famous Rayleigh-Ritz technique: 
• Restrict the operator H to the subspace Dd by orthogonal projection P onto Dd'- 



H = H 



:= PHP . 



• Let Ek, k = 0, 1, . . . , d — 1, denote all d eigenvalues of this restricted operator 22, 

H\xk) = E k \x k ) , k = 0,l,...,d-l , 

ordered according to 

Eq < Ei < . . . < Ed-i ■ 
Then the fcth eigenvalue Ek (counting multiplicity) of H satisfies the inequality 

k = 0,1,..., d- 1 . 



Eh < Ei. 



Assume now that this d-dimensional subspace Dd is spanned by some set of d linearly independent 
basis vectors k = 0,1, ... , d—1: Then the set of eigenvalues E may immediately be determined by 
diagonalizing the d x d matrix 

(<V<|%j)), i,j = o,i,. ..,d-i , 

i.e., as the d roots of the characteristic equation 

det((^:|ff|^)-^(^l^)) = 0, i,j = 0,l,...,d-l . 

(This becomes clear by expanding any eigenvector of the restricted operator H in terms of the basis 
vectors \ipk), k = 0, 1, . . . , d — 1, of the subspace Dd-) 



3.7 Variational Bound for the Ground State |T0| , ^] 

In order to derive, as a first application of the Rayleigh-Ritz variational technique, an upper bound on 
the ground-state energy eigenvalue Eq, it is sufficient to focus one's interest to the case k = (i.e., to 
consider an only one-dimensional subspace). In this particular case, the min-max principle reduces to 

(This inequality simply states that the ground-state energy Eq is less than or equal to any expectation 
value of the Hamiltonian H .) However, one is certainly entitled to consider simultaneously even sets of 
one-dimensional trial spaces and to compute optimized upper bounds by the following prescription 
(or "recipe"): 

1. Choose a suitable set of trial states {IV'a)} (with elements distinguished from each other by some 
variational parameter A). Each of these trial states serves to span a one-dimensional trial space. 

2. Calculate the expectation values of the Hamiltonian H with respect to these trial states 



25(A) 



(ipx\ip>< 



3. Determine (from the first derivative of E(X) with respect to A) that value A m i n of the variational 
parameter A which minimizes 25(A). 

4. Compute 25(A m i n ) (that is, the minimal expectation value of H in the Hilbert-space subsector of 
the chosen trial states { IV'a) }) - 

When going through these steps, your reward will be an optimized upper bound for the ground state: 

25q < 25(A m i n ) . 
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In order to get rid of the troublesome square-root operator in the Hamiltonian H , one may adopt a 
trivial (nevertheless fundamental) inequality for the expectation values of a self-adjoint but otherwise 
arbitrary operator O = O^ and its square, taken with respect to any arbitrary Hilbert-space state 
in the domain of this operator: 

(mow 2 < mo 2 m a>\i>) 

and thus 



mom ^ 


l(Mo 2 W 


- V 





It should prove to be advantageous to consider the above inequality for the kinetic-energy operator T: 



\mm < mm) 

Apply this to the semirelativistic Hamiltonian H: 

E < ^ (by the min-max principle) 

(ip\ip) 

W\T + V\iP) 

" VW^W (by the above mequahty) 



V MM MM ' 

Note, with due satisfaction, that now only the expectation values of p 2 and of V have to be evaluated. 

For the Coulomb potential (||), maybe the first choice for the coordinate-space representation of 
trial vectors which comes to one's mind are the well-known (normalized) hydrogen- like trial functions 

ipx(x) = \ — exp(— Xr) , A > . 

V 7T 

This choice of trial states yields for the expectation values 
• of the square of the momentum p: 

(i>\ \p 2 \iP\) 



of the inverse of the radial coordinate r: 

(^a|-0a) 



= A 2 



A . 



Insertion of these expectation values into the above boundary expression results in the one-parameter 
set of upper bounds 

E < VA 2 + m 2 - a X for all A > , 

with the absolute minimum 



Eq < m yl — a 2 



this is identical to the previously found generalized operator bound for n — 1 . For a^O this bound is 
lower and thus better than the Schrodinger bound on the ground-state energy (which is characterized 
by the quantum numbers n r = £ = and thus n = 1): 



E s .o = m 1 



a 2 



One may conclude that the variational technique yields indeed improved upper bounds on the energy 
levels. 
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3.8 Energy Levels at the Critical Coupling Constant fp7l| ] 

For the Coulomb potential (g), it might be of interest to derive upper bounds on the energy levels at 
the critical coupling constant a c , particularly, to improve the upper bound on the ground-state energy. 
A still rather simple-minded but very useful choice for the basis vectors \ijik) (labelled by some positive 
integer k, i.e., k — 0, 1, 2, . . .) of the d-dimensional trial space Dd is given 



in configuration-space representation by 



/ (2m) 2k+2 l 3+1 

^^ = V Jr(2fc + 2 /j + i) rWexp( - mr) ' r =W' m>0 ' 

in momentum-space representation by 



/ (9 777*12 fc+2/3+l Sill 



p 

(k + 8 + 1) arctan — 

m 



27T 2 r(2fc + 2/?+ 1) V H ' p(p 2 +TO 2)(fc+,3+l)/2 



/3 > 



Here, the parameter 8 allows, for a given value of the coupling constant a, of the total cancellation 
of the divergent contributions to the expectation values of, 

• on the one hand, the kinetic-energy operator T for large momenta p and, 

• on the other hand, the Coulomb interaction-potential operator Vc(r) at small distances r: 

8 = 3(a) . 

The parameter 8 is implicitly determined as a function of the coupling constant a, for instance, by the 
relatioij] g @ 



a — 8 cot y- fij f° r the ground state 



Hence, the critical coupling constant, a c , is approached in the limit of vanishing 8, that is, for 8 — > 0. 
For the present choice of trial states, the above-mentioned singularities arise only in matrix elements 
taken with respect to the ground state \ipo), i.e., in the matrix elements (ipo\T\ipo) and (V'ol^cMIV'o)- 
For arbitrary parameter value /3, the above choice of basis vectors \tpk) implies for the matrix elements 

• of the kinetic-energy operator T 

2 i+j+2p+i m T(i + 8 + l)T(j + 8+1) 



(V><|T|Vi> 



d?y 



T v /r ( 2 * + 2/3+l)r(2j + 2/5+l) 

cos[(i — j) arctan y] — cos[(i + j + 2 8 + 2) arctan y] 

(1 + y 2)(i+j+2/3+l)/2 





which has to be evaluated with the help of the expansion 

cos(iV arctan y)= {l + y2)N/2 E( n ) cos ) ^ for iV = 0, 1, 2, . 

of the Coulomb interaction-potential operator V = Vc (r) 

2mar(j + j + 2/3) 



v/r(2i + 2/3+l)r(2j + 2/3+l) 
and for the projections of the basis states \ipk) onto each other 

/AI/lN r(i+j + 2/3 + i) 



Vr(2i + 28 + i)r(2j + 2/3 + 1) 



4 Choosing equidistant values for the parameter 8, this relation then yields the curious values of the 
Coulomb coupling constant a adopted in Tables [l], |[ and [| 
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For notational simplicity, define a dimcnsionless energy eigenvalue e by extracting overall factors: 



E =: — me . 

7T 



The resulting characteristic equation for the case of relevance here, viz., /? = 0, is typically of the form 



det 



4 In 2 - 2 - e 


V2 


e 






3 


"71 






V2 e 


17 






= 


T~ V2 


15 


— £ 







V 



The roots of the characteristic equation (which may be calculated algebraically up to d = 4) then read, 
• for d = 1, 

£ = 2(2 In 2-1) , 
which entails, as upper bound for the ground state, 

— = 0.4918... ; 



for d = 2, 

£= 1~5 ( 601n2 ~ 23± ^(60 In 2) 2 - 4800 In 2 + 1649 
which entails, as upper bound for the ground state, 

— = 0.484288... ; 

m 

while, for d — 4, a lengthy expression for the roots implies, as upper bound for the ground state, 



E 

— = 0.4842564 . . 
m 



Interestingly, already for d = 2 our (analytical) bound lies well within the numerically obtained range 

0.4825 < — < 0.4842910 for a = a c . 
m 

3.9 Generalized Laguerre Basis 

Quite generally, upper bounds on eigenvalues may be improved by (suitably) modifying the trial space 
Dd, by enlarging it to higher dimensions d or by spanning it by a more sophisticated set of basis states. 
A — rather popular — choice of trial functions is based on the generalized Laguerre polynomials (x) 
(for the parameter 7): these are orthogonal polynomials 

• defined by the power series 



4 7) w = E(- 1 )* 



t=0 



k + 7 \ x 



k-t t\ 



and normalized, with the weight function x 1 exp(— x), according to 

00 

r( 7 + fc + i) 



J dxx^ exp(-x) 4 7) (a0 4^ 0*0 



k\ 



Skk' ■ 



In order to construct trial vectors of the subspace Dd, corresponding to states with orbital angular 
momentum t and its projection to, introduce two variational parameters, namely: 

• [i (with dimension of mass), and 

• (3 (dimcnsionless). 
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The chosen set of "Laguerre" trial states \ipk) is then defined by the configuration-space representation 

with the spherical harmonics 3^m(^) f° r angular momentum £ and its projection m, depending on the 
solid angle f2 and orthonormalized according to 

d&y£ m {ft>) yt'm'ity = Sw 5 mm i . 

The requirement of normalizability of the (Hilbert-space) trial states \ipk) imposes on the parameters 
(i and (3 the constraints 

fj. > , 0>~~- 

In this case, the configuration-space trial function ^^ m (x) satisfies the orthonormalization condition 

d 3 a;^ ita (x) il>k',t'm'{x) = s kk> Sm 8mm' ■ 
The momentum-space representation of these trial states \rpk) (obtained by Fourier transformation) is 



t , x / (2At)"+^+ifc! (-i)<|p|< A (-l)V + 

VmMPJ - I m n, /) , t ,, 1 ,I + i/ 2r /,,3\L 



r(2£ + 2/3 + fc + l) 2^+ 1 /2 r (^+|) ^ t! 



fc - 1 



r(2l + /? + f + 2)(2/j)' f 2£ + P + t + 2 + t 3 p 2 , v 

o ' ^~ >* + o> „2 I ..2 yemitlp) 



( p 2 +M 2)(2 Wt +2)/2 ^ 2 ' 2 ' " 2'p 2 +MV P; ' 

with the hypergeometric series F, defined by 

_ T{w) ^ T{u + n)T(v + n) 
{U,V ' ,W ' ,Z) ~ T(u)r(v) ^ T(w + n) n\ ' 

Then trivially, the momentum-space trial function tpk,em(p) satisfies the orthonormalization condition 

d 3 pTpk,im(p)Tpk',t'm'(p) = $kk> Sa> S mm > ■ 



3.10 Power-Law Potentials 113 



In order to be fairly general, consider an interaction potential of power-law form: 



V(r)=^2a n 



r bn , r = Ixl 



where the sets of (otherwise arbitrary) real constants a n and b n are only constrained by the necessary 
boundedness from below of the Hamiltonian: 

b n > — 1 if a n < . 

The matrix elements of this power-law potential V(|x|) are easily worked out: 



\l T(2£+2p + i + l)T(2£ + 2P + j + l) ^ (2/x) 6 ™ § § ( till 



x ( I + Jit 2/, )( 3+ ;i:^)r(2. +2 / J+t „ + t + ., + i). 

From this expression, the potential matrix V = ({ipi\V\ipj}) may easily be written down explicitly: if 

• considering, for instance, only radial excitations (£ = 0) and 

• choosing for the variational parameter f3 the value (3=1, 

one finds _ 

/ 3 -V3b n ... \ 



-V^&n 3 + b n +b 2 n ••• 

v ; : ••. j 
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3.11 Analytically Evaluable Special Cases |jT3 



Compared with the above calculation, an evaluation of the kinetic-energy matrix elements (i(}i\T\ipj) is 
certainly more delicate. However, there are several situations which allow for an analytical treatment. 

3.11.1 Orbital Excitations 

Restricting to the case i = j = but allowing for arbitrary values of the orbital angular momentum I 
means to consider pure orbital excitations. Merely for definiteness, fix, for the moment, the variational 
parameter (3 to the value (3 = 1. In this case, the matrix elements of the power-law potential operator 
simplify to 

while a rather straightforward calculation yields for the matrix elements of the kinetic-energy operator 

Aiming at analytical results, one should clearly try to get rid of the hypergeometric series F. There are 
several possibilities to do so: 



In the ultrarclativistic limit (realized for vanishing particle mass, m = 0), by use of the relation 

T(w) T(w — u — v) , „, 

F{u, v; w; 1 = ; \ ( for w + 0, -1, -2, . . . , K (w - u - v > , 

1 (w — U) l (w — V) 

the above kinetic-energy matrix element simplifies to 

2[T{£ + 2)} 2 



(Vo|T|Vo> = 



r(* + f)r(*. 

For instance, for the linear potential 

V(r) = ar, a > , 

minimizing the expectation value (ipo\H\ipo) with respect to the variational parameter fj, leads to 
the minimal upper bound 



In the limit of large orbital angular momenta (I — > oo), this optimized upper bound reduces to 

lim (mm(ipo\H\ip )) =%a(t+- A 

l^roa \H>0 J \ 4 

Consequently, it describes linear Regge trajectories]^] 

[E(£)} 2 cx I , 

which is in striking accordance with all other findings based on different considerations |l4|, [l^] . 

Another possibility — which will not be followed any longer — is to fix the variational parameter \x 
to the value n = m (which reduces the last argument of the hypergeometric series F to 0) and to 
take advantage of the relation 

F(u, v; w; 0) = 1 , 
in order to obtain the kinetic-energy matrix element 

(ore,),'';"'";' 1 ., 



5 The square of the exact optimized bound is already almost perfectly linear. Its deviation from the 
asymptotic linearity is monotonously decreasing for increasing angular momentum i. The maximum 
relative deviation (occuring clearly at the point I = 0) is 1 — 5 7r/16 = 1.8 %. 
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3.11.2 Radial Excitations 

Focusing one's interest exclusively to states with vanishing orbital angular momentum I (£ = 0) means 
to consider only radial excitations. The matrix elements of the kinetic-energy operator then typically 
take the form 



4/3+1 i- 2 ) t+s 

t=0 s=0 



r(2/3 + i+ l)r(2/3 + j + 1) 2tt ^ ^ ^ t\s\ 

r(0 + t + i)r(p + s + i)i ts 



1 + 2(3 
i - t 



J + 2(3 
3 - s 



with some integral Its which under certain circumstances may be calculated analytically. For instance, 
for fi — m, and for 2 (3 integer and thus (remembering the above normalizability constraint 2 (3 > —1) 
non-negative, i.e., 2/3 = 0, 1, 2, . . ., one easily finds 



2/3 + t + s + \t-s\ + l 



1-1 \ts\ 

E 

n=0 



2/3 + t + s + \t - s\ - n 



\t-s\ 
n 

/n 7r\ 
\~2~) 



T(2j3 + t + s + 



-1 2/3+t+s+2 

E 



cos 



2/3 + i + s + 2 
n 



- r(2/3 + t + S+ l--)cos(-) 



With this result, the kinetic- energy matrix T = ((ipi\T\ipj)) reads, for variational parameter (3 fixed to 
(3 = \, explicitly: 

/ 



T = 



64 

157T 



1 

^3 
7 



V3 
7 
11 

y 



V 



In this way, analytic expressions for the matrix elements {tpi\H \ipj) of the semirelativistic Hamiltonian 
H with power-law potential may be found. 

• Up to a trial-space dimension d = 4 these matrix elements are, at least in principle, algebraically 
accessible. 

• For larger trial-space dimension, i.e., d > 4, the energy matrix {{^i\H\ipj)) must be diagonalized 
numerically, without, however, the necessity to invoke (time-consuming) integration procedures. 

In order to estimate the quality of all the upper bounds derived, consider, for the following reason, 
the harmonic-oscillator potential 



V(r) =ior 2 



u > 



In momentum space the operator r 2 is represented by the Laplacian with respect to the momentum p, 



while the kinetic energy, nonlocal in configuration space, is represented by a multiplication operator. 
Consequently, the semirelativistic Hamiltonian H in momentum-space representation is equivalent to 
a nonrelativistic Hamiltonian with some effective (square-root) interaction potential. The equation of 
motion resulting in this way is then solvable with one of the numerous numerical procedures designed 
for the treatment of the nonrelativistic Schrodinger equation. For the comparison (Table |^) of these 
Laguerre bounds with results of a numerical solution of the ("momentum-space," in the above sense) 
Schrodinger equation, we employ some standard set of parameter values ]l6| : m = 2 GeV, u> — 2 GeV 3 
(as well as the choice /u = m, /3 = 1); we observe a rapid convergence of our upper bounds towards the 
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Table 4: Upper bounds (in units of GeV) on energy levels (labelled in usual spectroscopic notation) 
for a d x d energy matrix. (Analytically obtained results in italics.) 



State 


1 x 1 


2x2 


25 x 25 


Schrodinger 


IS 


12162 


3.9276 


3.8249 


3.8249 


2S 




8.1085 


5.7911 


5.7911 


3S 






7.4829 


7.4823 


4S 






9.0215 


9.0075 



4 SUMMARY 

The min-max principle and related simple theorems allow to derive rigorous statements on the energy 
levels predicted by a bound-state wave equation like the spinless Salpeter equation. In principle, these 
methods should work for all physical situations which may be formulated as an eigenvalue problem. If 
one is not satisfied with the outcome of this procedure, or if one is interested in the corresponding wave 
functions, one has to actually solve the equation of motion by some numerical approximation method. 
The most efficient among these is perhaps the "semianalytical matrix method" proposed in Ref. |l7j] . 
Nevertheless, there still remain questions of the significance |l(| of too naive relativistic improvements 
of some nonrelativistic potential models for the description of hadrons as bound states of quarks . 
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